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Abstract 

Due to the lack of simulation tools that take into account the actual geometry of complicated 
quantum Hall samples there are lots of experiments that are not yet fully understood. Already 
some years ago R. G. Mani recorded a shift of the Hall resistance transitions to lower magnetic 
fields in samples of a Hall bar with embedded anti-Hall bar by using partial gating. We use a 
Nonequilibrium Network Model (NNM) to simulate this geometry and find qualitative agreement. 
Fitting the simulated resistance curves to the experimental results we can not only determine the 
carrier concentration but also obtain an estimate of the screened gating potential and especially 
the amplitude and lengthscale of potential fiuctuations from charge inhomogenities which are not 
easily accessible by experiment. 
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The main application of the classical Hall effect lies in the simultaneous determination of 
the carrier concentration and the mobility. For small enough temperatures and pure samples 
quantum effects become important. In the quantum Hall regime experimental results for 
magnetic fields corresponding to plateau transitions depend on microscopic details like the 
potential landscape, which complicates interpretation. 

In experiments of the quantum Hall effect (QHE) a fixed voltage or current is applied 
to two metallic contacts, while additional contacts are used to measure longitudinal and 
transversal potential differences. For interpretation it is therefore essential to know the dis- 
tribution of the electrochemical potential in the sample. Genericproperties of this potential 
for simple geometries have been analysed in a series of papers jl-3|. The plateaus at mul- 
tiples of h/e^ of the transversal (Hall) resistance as a function of magnetic fie^ 
explained in terms of a simple modification of the Landauer model by Biittiker 
non-interacting electrons in an empirical confinement potential {g]. However, the resistance 
between plateaus depends on the detailed geometry, electrostatics of electrons or the poten- 
tial landscape generated by excess charges in the doped semiconductors in vicinity to the 
2-dimensional electron gas (2dEG). 

Electrostatics was investigated self-consistently in simple geometries 2|, |7| and shown to 
lead to the formation of alternating compressible and incompressible stripes. One concludes 
that in the plateau regime current in response to (transversal) electric field only flows in 
the incompressible stripes. The picture becomes more complicated for nonideal contacts 
due to potential barriers or when applying gatings because channels are (partially) blocked 
and the electrochemical potential is changed in their surroundings. F. Dahlem et. al. found 
experimentally that width and magnetic field values of the transition region between plateaus 
can change significantly 18|. At the same time we investigated such situations with the 

fin 

nonequilibrium network model (NNM) [9|, 110[, with good agreement to the experiments. 

In addition the NNM has proven successful for exotic sample geometries like anti-Hall 
)ar s within Hall bars when compared to experiments of ungated samples by R. G. Mani 
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13| . Shortly after, Mani applied partial gating to his samples and recorded a shift of the 
Hall resistance transitions to lower magnetic fields as a function of the gating voltage fl^ . 
In the present paper we describe simulations of these samples with the NNM and show that 
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the features in the experiment are well reproduced. 

Moreover, by manually fitting the transversal resistance we extract the electron concen- 
tration, the screened gating potential and the average curvature of saddle points of the 
potential. From the curvature and a statistical model we obtain estimates for the amplitude 
and lengthscale of potential fluctuations. 



II. THEORY 



The exact Hamiltonian of a typical integer quantum Hall sample is complicated due an 
interface of two semiconducting layers, excess charges, and the presence of metallic con- 
tacts. Accurate calculations for the whole structure by ab initio methods fail due to the 
prohibitively high number of electrons in the heterostructure. 

For sufficiently high magnetic field electron wavefunctions are highly confined to the mag- 



netic length Ib = ^/h/eB, where e and B denote elementary charge and magnetic induction 
normal to the 2dEG. Therefore we model the potential landscape in the 2dEG phenomeno- 
logically, using a network of semiclassical trajectories in a (slowly varying) random potential 
that mimics the effects of the excess charges. The most prominent of this type of models is 
the Chalker-Coddington model 15|, which is able to predict statistics of states and scaling 
exponents but cannot describe the nonequilibrium steady state. In contrast, our NNM is 
designed to calculate these nonequilibrium quantities. The model rests on the local equi- 
librium approximation [16], which is applied to a network of semiclassical wavefunctions. 
In this way we attribute unique thermodynamical quantities such as the electrochemical 
potential to each single wavefunction. 

Transport by highly localized quantum states can be viewed as a percolation problem, 
where it is well known that only pivotal edges are relevant to the (bond) percolation problem 



171]. These edges correspond to saddle points of the potential landscape, which then form 
the network [15] . At each saddle 4 chemical potentials meet and we assume that phases are 
destroyed by decoherence, such that the distances between saddles become irrelevant. In 
this approximation we can replace the potential fluctuation by the model potential V{x, y) = 
V[cos{ijjy) — cos(co'x)] of a regular grid of saddle points, where the period L := 2h/ijJ and 
amplitude V should be understood as averages of the microscopic potential. 



3 



We define the ratio of longitudinal to transversal field component at a saddle point as 



„ E-r 
p ._ ±_ 

Ey Ui - U4 U2- Ms 

In this way we construct a "transfer" equation for the chemical potentials 
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The chemical potential distribution can be calculated as a boundary value problem once the 
values for P at each node are given. We model external contacts supplying current in the 
NNM by saddles with a pair of trajectories that point into/out of the sample. One of the 
two states is fixed to the chemical potential of the current supply contact. 

We neglect nonlinear effects, that is, the values of P do not depend on the chemical 
potentials. This should be well justified in case of typical experimental currents of yuA, 
corresponding to voltage differences lower than 0.1m V across the sample or energies well 
below the Landau level (LL) spacing of heB/m* ^ 1.72'imeV at B = IT, respectively. 

Except in a small region around the hotspots (at opposite edges of the current inducing 
contacts) the local conductance tensor can be approximated as purely off-diagonal. This 
leads immediately to P = ly/Ix- This approximation has the advantage that only the electric 
field has to be calculated. According to the edge channel picture we call T the probability 
of transmission in longitudinal direction, therefore we get ly (x R and cc T = 1 — R. P 
can then be calculated from elastic tunneling transition probabilities at the Fermi energy 



across a saddle as 
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P = ^mn — = exp -e— , (3) 

with e := Ep — En — Vs the energy of the trajectory, in terms of the difference of the Fermi 
energy to the energy of state n, relative to the saddle energy Vs- Vs reflects the sum of all 
potentials (disorder, confinement...). We define the "center of the transition" as e = 0. While 
B denotes the magnetic field strength, the parameter c := ^ is related to the curvature 
of the model potential energy fluctuations with amplitude V. Similar to other approaches 
[sl, [2^ we calculate longitudinal and transversal resistance from the potential distribution 
by identifying the current with the macroscopic current direction. 

We use LLs and a simple self-consistent Thomas-Fermi approximation where the po- 
tential at zero temperature is obtained from the charge density by multiplication with 



C = 50mV^/(10^^cm~^) (units e = 1). Furthermore we use a constant broadening of 0.5meV 
to mimic scattering by a short range disorder potential. Equihbration among edge channels 
is taken into account by assuming tunneling and using an exponential function with a decay 
parameter. 



III. RESULTS 



A. Transversal resistance and fit of experiment 

Figure [1] schematically shows the geometric setup of an anti-Hall bar embedded in a Hall 
bar, used in the experiment [14]. The gated rectangular region is filled with yellow color and 
a confinement potential develops near the inner and outer border of the 2DEG, indicated 
by solid lines. In the experiment a constant inner/outer current was driven through the 
device by applying appropriate voltage differences at the contact pairs A-B of the outer 
Hall bar and 1-2 of the inner anti-Hall bar. We denote in the following the longitudinal 
(parallel to direction A-B) and transversal direction by x and y, respectively. Transversal 
Hall resistances Rxy are obtained from F-C and 6-3. As discussed in [ol it is sufficient 
to use point contacts in simulations of macroscopic Hall samples, such as in the present 
experiments. 

We manually fit our simulations to the experimental results in order to obtain important 
microscopic information like the carrier concentration, the enhancement factor for the Zee- 
man energy, the screened gating potential, and the magnitude and correlation length of the 
potential fluctuations. We estimate errors from half the thickness of lines in plots of the 
experimental results. 

To get the carrier concentration we fit the plateau transition center of R^y of the outer 
ungated Hall system for the (3 last integer) transitions from filfactor z/ = 4toz/ = 3, z/ = 3 
to z/ = 2, and z/ = 2 to z/ = 1. This results in a carrier concentration oi n = 1.81 * 10^^ ± 
0.02 * lO^^cm"^. Our value is slightly below the range 2 * lO^^cm"^ < n < 3 * lO^^cm"^ 



2l| . where however a fit to the classical Hall slope is 



proposed in the experimental work 
normally used. 

Due to spin polarization a Zeeman term adds to the energy, resulting in an energy 
difference of AEz = gfiB/2, with g an enhancement factor and := gcaAsfJ'B- Here 
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2.20 


2.26 



TABLE I: Magnetic fields in T of the transition centers for various bare gating potentials. We 
estimate the error as 0.05T. 

QGaAs = —0.44 is the Lande factor of GaAs heterostructures and fiB = eh/ 2m* the Bohr 
magneton for electrons with effective mass m*, which we set to the typical effective mass 
of GaAs, m* = 0.067m. By manually fitting the transition centers in case of no gating we 
arrive ai g = 12.5 ± 1.5. Such large enhancement values are typically seen in transport 



measurements 
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We fit the magnetic field at the center of the plateau transitions observed in the inner 
and outer leg of the anti-Hall structure, and arrive at bare gating potentials of lOmeV, 
SOmeV and 50meV^. The screened gating potentials Vg{B) as functions of the magnetic field 
B are shown in Fig. [3] for these bare gating potentials. We note that within the plateau 
transition Vg is small and fairly constant only for small values of Vg. The magnetic fields at 
the transition centers are then collected in the table [H 

Figure E] shows Rxy for various Vg applied to the inner Hall bar. The dominant feature 
of gating lies in a shift of the plateau transitions to lower magnetic fields due to a reduced 
electron concentration. We note that the Fermi energy is fixed by a reservoir and therefore 
only a function of the magnetic field but not the local gating. The impression of a more 
or less rigid shift is a consequence of Ef{B) — Vg{B) having small variations on the scale 
of the magnetic field interval of Rxy transitions. We mention that if the gating is large 
enough, which is the case for V^'' = 50meV^, the center can jump down to the next lower 
linear increasing part of the Fermi energy, corresponding to the LL band at lower energy. 

In order to obtain the potential curvatures for all gate potentials, we would have to fit the 
curvature of each transition and for each Vg separately, which is very demanding. Therefore 
we determine cq using the ungated Rxy curve and calculate from it the respective value Cy for 
Vg. We focus on the three transitions with lowest filfactor u, starting with center at highest 
magnetic field and ordered due to decreasing magnetic field (corresponding to the states 
(i^, s) G {(2, +), (2, — ), (1, +)} crossing the Fermi energy, where s is the spin orientation). 
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TABLE II: Fitted curvatures cq in units meV/T for various transitions and bare gating potentials. 
We estimate the fitting error as ziz0.2meV * T from the hne thickness of the Hall resistance in the 
experimental figure. 
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TABLE III: Curvatures cq in units meV/T for various transitions and bare gating potentials for 
the part of the anti-Hall bar with no gating applied 

We summarize the manually fitted cq values as a function of Vg in table Ull 

It is interesting to also present the cq values at the leg without gating, summarized in 
table mil For each transition the values should not depend on Vg, if there is no charge 
transfer from the gated to the ungated leg. It seems that the curvature decreases (the 
potential landscape is getting flatter) with increased gating on the opposite leg. On the 
other hand, we expect the charge transfer to be small in such a macroscopic bar. Within the 
(large) uncertainties of the experimental curves we cannot predict even a qualitative trend. 
Therefore we average along each row to obtain 0.17, 0.4 and 0.67 in units of meV/T for 
4 — 7- 3, 3 — 7- 2, and 2 — 1, respectively. 

We stress that the shift of the transition center with gating also changes the slope of 
Rxy in the center of transitions (e = 0) due to the appearence of B/c (note c = a/2, to be 
consistent with the definition of a in |9|]) in P (see Eq. [3]). In order to compare transitions 
for various gatings, we demand that the interval in magnetic field 6 By, corresponding to 
\e{B)B/ci,\ = 1, is the same as the one at zero gating, given by \e{B)B/co\ = 1, using 
an appropriate definition of Ct,(co). In order to get explicit expressions we expand linearly 
around the magnetic field at transition centers, e{B) ~ eo + eiBtr- In this way we get from 
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TABLE IV: Effective curvatures Cy in units meV/T for various transitions and bare gating poten- 
tials used in the experiment 



\e{B)B/c\ = 1 the interval 




(4) 



Demanding 6By = 6Bq we map the two transitions onto one curve and arrive at the curvature 
of the gated sample, 
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Using the values of Cq in table [TTl we calculate for each transition in terms of Eq. |5l 
which we summarize in table IIVI 



B. Potential landscape 



To determine both the amplitude and lengthscale of the potential fluctuations we need 
another quantity besides the curvatures. In this resp ect a model of charge fluctuations by 



As is well known, if a collection 



statistically independent electrons is very useful [23 
of statistical objects with identical properties (described by the same random variable) are 
independent then the relative fluctuations of X := Xj are given by AX/ < X >oc N"^^"^, 
where < X > denotes the average of X and AX : = 
XjXk > - < X >2]VV^ together with < X^X^ >= 6jk 



'< (X- < X >)2 > 



E". < 



We partition the sample in N cells of equal size and interpret the electron distribution (or 
number of electrons) in each cell as a random variable. Assuming charge neutrality on the 
average the absolute fluctuations of the charge density n are then 5n = n/N^/"^. Employing 



the Poisson equation we then get the estimate VhneV] 



the (long range) potential energy fluctuations |23|. n = 1.3797 * 10 ^C/Vm denotes the 
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for the magnitude of 
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TABLE V: Number of correlated electrons in statistically independent cells. We rounded to the 
next higher /lower integer. 

dielectric constant of the sample (GaAs) and K is the modulus of the smallest wavevector 
supported by geometry. 



The detailed description of the experimental setup [IJ] lets us estimate the length of 
each leg as / = 2mm and the transversal width as w = 0.2mm. Assuming variations of 
charge only normal to equipotential lines the largest wavelength should occur in the middle 
of the plateau transition where the Hall angle is close to it /A. We arrive at a wavelength 



of 2tc/K = 2 * 10~^\/2 meters. Moreover we get 6n = n/y/N = ^JnNc/lw. Noting that 
n hardly varies between transitions we use averages for each V^. The fitted densities (in 
lO^^^cm^^) turn out to be ri = 1.81, n = 1.62, n = 1.25, and n = 0.88 for bare gating 
(in nieV) Vg = 0, V^" = 10, Vg = 30, and Vg = 50, respectively. Nc denotes the number 
of (correlated) electrons in each fictitious cell of the statistical model. The amplitudes of 
potential fluctuations become V j \fWc = 2.203, 1.972, 1.521, and 1.071, respectively in units 
meV, for increasing gating. 

This gives the interesting possibility to get a value for the number of correlated electrons, 
as explained in the following: A plausible upper bound for is given by the energy difference 
between LLs^because the screening tries to suppress higher potential amplitudes due to 
Wulf et. al. [23|. Namely, if the potential energy exceeds the LL spacing a new LL band 
is occupied, which leads to strong screening that suppresses the potential until the new 
LL band is emptied and screening is weak again. Clearly this argument does not hold 
for potentials that are so large as to change the level structure significantly or lead to the 
breakdown regime. We arrive in this way at the upper bounds A*"" that are presented in 
table |Vl The variation with Vg turns out to be small. 

Using the definition of the curvature we get the potential correlation length in nm as 



transition 



2^1 138.8 ± 55.5 138.0 ±27.1 125.3 ± 13.5 164.9 ± 24.6 
3^2 148.2 ±87.1 (190.3 ± 111.9) 148.2 ±31.5 179.5 ± 44.c 
4^3 (166.5 ± 151.3) 160.1 ±66.7 227.1 ± 113.6 (275.3 ± 162.0) 

TABLE VI: Correlation length and their errors of potential energy fluctuations in nm, using the 
values of Nc from table |Vl Brackets indicate values with large errors. 



where V 



ev 



^Nc/h with numerical values v = 95.46, 85.44, 65.93, and 46.41 for increasing 



gating. The so obtained potential correlation lengths are summarized in table IVll We judge 
the error as coming predominantly from the curvature, because the concentration n can be 
fitted with high accuracy. The resulting error varies strongly due its proportionality to c~^/^. 
Qualitatively the error is less for less gating and smaller fillfactor (that is, for larger value 
of L inlYlD. 

Table Wl\ which shows microscopic information that is hard to be addressed by experi- 
ment, is our main result. These values should be viewed as the largest correlation lengths 
appearing in the potential landscape, due to the Hartree-like estimation of the potential am- 
plitude used. The general trend is decreasing screening (smaller L) with increasing Vg. This 
makes sense, despite the high error for low L, and is due to decreased carrier concentration. 
Similar the screening increases with the number of electrons present at transitions, that is 
with decreasing magnetic field, due to the same reason. However, there are exceptions to the 
rule which seem to come from the non-trivial dependence of the Fermi energy and screened 
gating potential energy with the magnetic field. 

We stress that via K a larger area of the 2dEG sample leads in principle to higher 
potential fluctuations, which gives different A^^^ due to the saturation at the LL energy 
difference. Furthermore, the presented (large scale) correlation lengths have a realistic order 
of magnitude when compared to experiments 25l-l27l| as well as to our own Hartree-Fock 
(HF) calculations, using a refined version, with Broyden mixing for fast convergence, of a 



code originally written by Romer et. al. [28|. We note here that in reality due to imperfect 
screening a reminescence of the distribution of charge from doping centers near the 2dEG 
affects the assumed statistical independence slightly. However, as shown by Gudmundsson 



et. al. 2J] the agreement with experiments is good. 
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C. Chemical Potential distribution 

Finally, in Figure [5] we present details of the electrochemical potential distribution for 
the case corresponding to a bare gating of Vq = SOmeV (— 300m\^ in the experiment). 
We selected 4 magnetic field values to show the generic behaviour. We first note that with 
increasing B the maximum voltage difference in the system increases, which is a consequence 
of constant injected current. At 3.66T the inner (gated) bar is in the transition region, which 
can clearly be seen from the longitudinal gradient in blue/red color in the upper /lower leg. 
On the other hand the red color on the other bar is homogeneous, that is, in the plateau 
and dissipation only occurs at the current contacts. At 4.20T inner and outer bar are in the 
plateau, therefore only a transversal gradient develops and the line of zero potential (white) 
is parallel to the longitudinal direction along the leg. 4.86T shows qualitatively the same 
picture as 3.66T with the role of inner and other bar reversed. Now the outer bar shows a 
longitudinal gradient. Finally at 5. SOT both bars are in the plateau again. 



IV. SUMMARY AND CONCLUSION 



In summary we present simulations of experiments done on an anti-Hall bar within a Hall 
bar geometry by Mani jl^ , which show a shift of the Hall resistance to lower magnetic field 
by applying partial gating. We calculate the chemical potential distribution for the integer 
quantum Hall regime under a constant current condition, as used in the experiment. The 
so obtained transversal resistances as a function of magnetic field compare well with the 
experimental curves. 

Fitting by hand the position of plateau transitions of the transversal resistance for applied 
gating, we are able to obtain the electron density, the enhancement factor for the Zeeman 
interaction, and the screened gating potential in the 2d electron gas. Finally we determine 
the curvature of potential saddles by fitting the width of plateau transions using the ungated 
Hall resistance curve. Employing a model of partitioning the sample into statistically inde- 
pendent cells with correlated electrons in each cell, we arrive at the amplitude of potential 
fluctuations. Defining a transformation we obtain plateau transition intervals for the gated 
system from the fitted values of the ungated system. In this way we are able to provide 
the magnitude and lengthscale of potential fluctuations from charge inhomogenities as a 
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function of fillfactor of the transition and bare gating potential. 
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FIG. 1: (color online) Geometry of the anti-Hall bar within a Hall Bar: The figure shows the 
labelling of contacts used and the shaded, yellow area denotes the partial gating applied as in the 
experiment. 
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FIG. 2: (color online) Variation of the transversal resistances as a function of magnetic field: The 
upper plot shows Rxy for various gating potentials applied to the inner gated bar. The lower curves 
show transversal resistance curves for various curvatures c and no gating. 
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FIG. 3: (color online) Screened gating potentials as a function of magnetic field: We show results 
for bare gatings of lOmeV, SOmeV, and 50meV. Low/high screening does not directly correspond 
to transition/plateau regions of Rxy as broadened LLs for different spin significantly overlap for 
the Zeeman energy found. 
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FIG. 4: (color online) Fermi energy and 4 Landau levels of lowest energy (2 spin-resolved pairs) as 
functions of the magnetic field. 
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FIG. 5: (color online) Distribution of the electrochemical potential: we show results for a bare 
gating potential of SOmeV and magnetic field values (from top, left to bottom, right) of 3.66T, 
4.20r, 4.86T and 5. SOT. 
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